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Gluon contributions to the pion mass and light cone momentum fraction 
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We calculate the matrix elements of the gluonic contributions to the energy-momentum tensor 
for a pion of mass 600 < M„ < 1100 MeV in quenched lattice QCD. We find that gluons contribute 
(37 ± 8 ± 12)% of the pion's light cone momentum. The bare matrix elements corresponding to the 
trace anomaly contribution to the pion mass are also obtained. The discretizations of the energy- 
momentum tensor we use have other promising applications, ranging from calculating the origin of 
hadron spin to QCD thermodynamics. 
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Introduction. — A striking feature of QCD is the large 
contribution of gluons to the mass and momentum of 
hadrons, so it is of fundamental interest to calculate the 
contributions of gluons from first principles using lattice 
QCD. 

The first moments 



<*>f(9 2 ) = 
<z)g(9 2 ) = 



£/=M, s ftxdx {/Or, g 2 ) + f(x, q 2 )} (1) 
Jgxdx g{x,q 2 ) (2) 



of the quark and gluon distribution functions fix), f{x) 
if = u,d, s, . . . ) and <?(x) acquire a precise field-theoretic 
meaning via the operator product expansion in QCD. 
They satisfy the well-known momentum sum rule (MSR) 
(x)f(g 2 ) -I- {x) g iq 2 ) = 1 and are related to the correspon- 
ding contributions to the energy-momentum tensor T^ v 
evaluated on the hadronic state. Separating the traceless 
part T pv from the trace part S for gluons, denoted 'g', 
and quarks, denoted T, T pv has the explicit form 



Tfj, u = T + T + i<5 M „(5 s + S l ) 



T s — i(5 F a F a 
fiv ~ i u t^ r pa r pa 

T l = 



4^ 

pa pa 
fia. vol ' 



(3) 
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where D p =D p — D^, Pig) is the beta-function, 7 m (g) is 
the anomalous dimension of the mass operator, and all 
expressions are written in Euclidean space. For an on- 
shell particle with four-momentum p = (i_E p ,p), = 
M 2 + p 2 , we have the relations 



(*, P |/d 3 zr; s (z)|vi/,p) 

(*,p|/d 3 z^(z)|vl/,p) 

(x) { + (x) g 



[E p - \M 2 /E p ] (x) t , g , (4) 
iM 2 /E p ) bt, s , (5) 
bt + b g = 1, (6) 



where states are normalized according to (p|p) = 1. We 
shall return to the renormalization of (a;)f, g below. 

Equation|4]shows that in the infinite momentum frame, 
where E p ~ P — > oo, (x) g represents the momentum 
fraction arising from gluons, and calculating (x) g is the 
main goal of this work. In the rest frame, the gluon 
contribution of Eq. [4]to the hadron mass is jM(x) g 
From Eq. [5] in the rest frame, the contribution of the 
trace anomaly S 9 to the hadron mass is \b g M and 
in this work we perform the first step to calculate this 
matrix element as well. 

Whereas non-singlet matrix elements can now be cal- 
culated to high precision in full QCD in the chiral 
regime @, [U, 0] , calculations of matrix elements of singlet 
operators are far less developed due to the computational 
challenges of calculating disconnected diagrams, which 
require all-to-all propagators, and matrix elements of 
gluon fields, which are notoriously noisy due to quantum 
fluctuations. The first attempt to calculate the quark mo- 
mentum fraction was in the proton in [5| , and was found 
to be numerically very challenging. In this exploratory 
study we treat the case of "heavy pions" with masses in 
the range 600 MeV < M n < 1060 MeV, where hadronic 
matrix elements in the quenched approximation, which 
neglects quark loops, are generally close to those in full 
QCD. The techniques developed here are applicable in 
full QCD calculations, and to the case of the proton. 

Lattice formulation. — We use the Wilson gluon action 
^E I , fl ^ Tr '{! - P^ix)}, where P pv is the plaquette, 
and the Wilson fermion action [|| at an inverse cou- 
pling 6/g 2 . = P — 6.0, corresponding to a lattice spacing 
a = 0.093 fm for ro = 0.5 fm There are two distinct 
ways Q to discretize the Euclidean gluon energy opera- 
tor T„ = 5(-E Q • E a + B a ■ B a ) and the trace anomaly 
S g _ ^M(E a • E a + B° • B a ) on a hypercubic lattice. 



'Electronic address: meyerh@mit.edu, negele@lns.mit.edu 



The first, denoted 'bp' for bare-plaquette, uses a sum of 
bare plaquettes P pv around a body-centered point Xq = 
x + i<i fi, which, when summed over a time slice, 



2 



yields 



2 X bp (<?o)Z g (5o) 



E 



(7) 



RcTr [E fc i%*(») - E k<l \\Pki{*) + Pki{x + aO)] 
a 3 yS^(x Q ) = 2x ' P(go) ff°~ 2 V ReTr x 

X X 

[E fe (l - Pok(x)) + £ fc<i (l - \[Pu{x) + Pki(x + aO)]) 
The other form, denoted 'be' for bare clover, is 

T* >)=* bC(go) f s(go) ReTr 



■9o 2 



5 bc (x)^ Xs bc (5o) 



dloga 



ReTr 



J2k( F 0k) 2 ~Y,k<l( F kl) 



where F liv {x) is the clover-shaped discretization of the 
field-strength tensor (see This form allows for the 



discretizations of off-diagonal elements of T M „ as well. 
Each of the normalization factors Z g (g ), x bc (ffo) an d 
X bc (5o) in Eq. d70) is of the form 1 + 0(g%). 

An additional freedom in discretization is local 
smoothing of the fields by replacing each link in Eqs. (|7l8p 
by a sum of a connected product of links joining the same 
two lattice points. This only changes the fields by higher 
dimension operators, and HYP smearing [11] is particu- 
larly suited for this application because it preserves the 
symmetry between all Euclidean directions and is lo- 
calized within a single hypercube. We use the original 
HYP-smearing parameters [111, and project onto SU(3) 
as in [l2l ]. 

Our criteria for the choice of the discretization are to 
maximize the signal-to-noise ratio, minimize cutoff ef- 
fects, and preserve locality as much as possible. The 
noisiest quantity we calculate is Tqo(x), which involves 
the near cancellation of E 2 and B 2 . Hence, we stud- 
ied the signal-to-noise ratio for four different discretiza- 
tions by comparing the variance of a related thermo- 
dynamic variable, the entropy density at temperature 
T = 1/Lq = 1.21T C Q, which is proportional to the ex- 
pectation value of J2 X Tqo(x), on an L x L 3 lattice with 
L/a = 16 and Lq/cl — 6. The resulting variances for the 
plaquette and clover discretizations with bare and HYP 
links are shown in Table |TJ We find dramatic differences 
between the discretizations, with HYP smearing reduc- 
ing the bare plaquette variance by a factor of 41 and the 
HYP-clover operator reducing the variance by a factor 
of 87. Variance reduction comes at the cost of a certain 
loss of locality, since the HYP plaquette and HYP-clover 
operators have extent 3a and 4a respectively. 

The normalization factor Z g (go) appearing in Eq. [7] 
is dictated by an exact lattice sum-rule for the Wilson 
gauge action and is known with a precision of about 1% 
(see [l3[ and Refs. therein). To obtain the absolute nor- 
malization of other discretizations, it is sufficient to com- 
pute their normalization x(do) relative to that of the bare 





relative variance 


normalization 


bare HYP 


bare HYP 


Too 


plaq. 
clover 


26.4(71) 0.6518(43) 
3.85(11) 0.3049(41) 


1 0.5489(68) 
2.184(67) 0.613(20) 


S 


plaq. 
clover 


2.64 (12) 0.474(13) 
1.180(39) 0.2975(72) 


1 0.9951(77) 
4.062(30) 1.410(13) 



TABLE I: Left: the relative variance, (£?)/{0) 2 ~ 1, of the 
operators O = ^2 x (o(x) — (°)°) (top: o = Too, bottom: o = S) 
on a 6 x 16 3 lattice at f3 = 6.0 for different discretizations 



described in the text, 
(top) and Xs{go,a/L ] 



Right: the normalization x(go,a/Lo) 



(bottom) of the operator relative to 
the bare plaquette, determined on the same lattice. 



Wilson gauge action P=6.0, (L /a) xl6 3 lattice 
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FIG. 1: A study of cutoff effects: the normalization 
x(go, o,/Lo) of three discretizations of Too relative to the one 
based on the bare plaquette as a function of Lo/a. 



plaquette, and the resulting \' s are given in Tab. U for 
the four discretizations. 

As a compromise between locality and variance reduc- 
tion, from now on we work with the HYP-plaquette dis- 
cretization. We performed a check of its discretization 
errors by computing the dependence of \ on a/Lo; which 
is a nonlocality effect. Figure [1] shows that the depen- 
dence of x on a/Lo is mild and statistically consistent 
with zero for L /a > 6, and that all four lattice oper- 
ators are viable discretizations of the same continuum 
operator. As a check of the correct normalization of the 
chosen HYP-plaquette operator, we computed its expec- 
tation value on the lightest scalar glueball. In that case, 
we know that the momentum fraction carried by the glue 
is one (see [l4| for an early calculation in SU(2) gauge 

theory), and indeed we find (x))r' = 1.16(18). 
The gluon momentum fraction in the pion. — We con- 
sider a triplet of Wilson quarks, labeled u, d, s, with peri- 
odic boundary conditions in all directions and with com- 
mon k = 0.1515, 0.1530 and 0.1550 corresponding to 
pion masses approximately 1060, 890 and 620 MeV on 
lattices 32 • 12 3 , 32 • 16 3 , 48 • 16 3 and 24 4 . To calculate 
the gluonic momentum fraction in the pion, we define the 
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Wilson action p=6.0 K=0.1515 



<x> g>efl (x ) 



32xl2 J 
32xl6 3 
48x1 6 3 
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FIG. 2: The effective gluonic momentum fraction, Eq.|9] in a 
heavy pion, M n ~ 1060MeV. 



effective momentum fraction 
8 a 3 



<*>2 ff w in ) 



3M W |A | 



(9) 



x; x eA 



E y Q-(Q) roo P (^)j(^,y)) ^ hp . 
E y ^(Q)i(^,y')> ( oo(a;0)) 



and similarly for 6 g bare ' ) by substituting T q 

min L Q mjn_, ) £o i ^min 

Lq , . . . , 2 j^q u, 2 -f Xq , . . . 



S hp . Here 

A _ / min Lo ^,min Lo i ^.min r min _ 

JlO — t-^o > • • • ) IT u > "2 ' > • • • > n x 

a}. This corresponds to creating a pion at the origin, 
annihilating it at the middle time slice, measuring the 
gluon operator over all times at least x nin away from the 
source or sink, dividing by the corresponding pion two- 
point function, and subtracting the vacuum expectation 
value of the operator. For large L and a;™ 111 , — > 

As a source field for the pion, we use the isovector 
pseudoscalar density j(x) = d(x)"f^u(x). Its two-point 
function is positive on every configuration, for each of 
which we do 12 inversions corresponding to Dirac and 
color indices. On a 24 4 lattice, we take advantage of the 
symmetry between all directions to perform these inver- 
sions at the points k(6, 6, 6, 6) for k = 0, 1, 2, 3 and sym- 
metrize expression ^ with respect to all directions, so 
that J2 X u^mm( x ) vanishes on every configuration. Fig- 



ure [2] shows our stable plateaus for (x)^ J ff at large values 
of x™ m for two lattice sizes, and all the results are sum- 
marized in Tab. HH 

Equation [5] has been derived for QCD at finite lattice 
spacing in [15| . In particular we have 



1 = (x) g + (x) f , <*) f = Z f (g )(x)^ e , 



(10) 



where, disregarding disconnected diagrams, (x} b arc has 
been computed in [201 ] at the same bare parameters (/? = 
6, k = 0.1530). The factor Zf(go) is the fermion analog 
o{ Z g (g ), see Eq. Q. 



Mir (Mev) 


32 • 12 3 


32 • 16 3 


48 ■ 16 3 


24 4 


1060(10) 
891(9) 
624(6) 


0.39(6)23091 


0.29(9)7ii3 


0.40(8)8331 


0.34(9) 1048 

0.36(8)3066 
0.58(16)2538 


1060(10) 
891(9) 
624(6) 


0.89(3)23091 


0.95(5)7113 


1.00(4) 8 331 


0.77(7)io48 

1.02(6)3066 
2-2(1)2538 



TABLE II: The glue momentum fraction (x)^ (top) and the 
bare trace anomaly matrix element bg barc ' (bottom) in the 
pion. The integer in each subscript denotes the number of 
configurations used. 



Renormalization of (x) g . — Recall that, in QCD, the 
renormalization pattern in the singlet sector reads [l6| 







_ t oo(m) _ 





Z sg 1 - Z s 
1 - Z„ g Z B 



r oo(ffo) 
T oo (So) 



(11) 



provided T 00 (go) are normalized so that Eqs. (|4I6[) 
hold. In lattice regularization, this requires the scheme- 
independent Z g (go) and Zf (go) factors, while Z gg and Zg 
are scheme-dependent functions of (afj,,go). The renor- 
malization group equation then takes the form 



fldu 



(z)g(M 2 ) 



gg(s)_ 

- c gg(sO 



-Cff(ff) 

Cff(ff) 



(^)g(M 2 ) 
<*}f(M 2 ) 



1] = -9 2 [£m 
respectively [13, H 
the linear combination 



+ c ff ] and c gg: s(g = 
. Besides the zero- 
[1 + rifxWlM + 



with fidfi log[Z gg + Z s 

^> 12tt 2 ' 9tt 

mode Too 

t((i)T 'oo(a0 renormalizes multiplicatively with anomalous 
dimension — g 2 [cg + c gg ], where /id^T = —g 2 [(cs + c gg )T + 
Cfj] . Note that the asymptotic glue momentum fraction is 
given by c ff (0)/[c ff (0) + c gg (0)] = ^ gg (oo) = 1-Z ff (oo) = 
-t(oo) = 16/[16 + 3iVf]. 

In the quenched approximation, Z gg — 1 due to the 
absence of quark loops [U [53]. This implies that the 
singlet part renormalizes multiplicatively and with the 
same anomalous dimension as the non-sing let p art, which 
has been computed non-perturbatively in [l9| . 



(^)g(M 2 ) 

<A(m 2 ) 



(x) s + [l-Z s (a(i, 9o )] (x} { , (12) 
Zs(a/J,,g ) (x)[. (13) 



The factor Zf(go) = 1 + O(g^) is as yet unknown beyond 
tree level. If we allow for a conservative error, based on 
the typical size of one-loop corrections, Zf(go) = 1-0(2), 
then using (x} baro = 0.616(4) and Z s (afi, g )Z f (g ) = 
0.99(4) for the MS-scheme at fi = 2GeV 0, [lj, our 
final result is 

( x )W(p2_ = 4GeV 2 ) = 0.37(8)(12) (M w = 890MeV), 

where the first error is statistical and the second comes 
from the uncertainty in Zf(go). 
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Finally, our result for the glue momentum fraction in 
a (heavy) pion is compatible with phenomenological de- 
terminations m>m> = °- 38 ( 5 ) at q 2 = 4Gev2 : 

based on Drell-Yan, prompt photoproduction, and the 
model assumption that sea quarks carry 10-20% of the 
momentum. The agreement suggests a mild quark-mass 
dependence, but only a calculation in full QCD and 
at smaller masses can substantiate this. Our result at 
Q 2 = 4GeV 2 lies clearly below the N{ = 3 asymptotic 
glue momentum fraction of 0.64. The fact that our result 
and the valence quark momentum fraction, computed 
in m^ a dd up to 0.99(8)(12) suggests that the omitted 
disconnected diagrams are small. 

Discussion of b g . — In a chirally symmetric formulation 
of massless QCD, the trace anomaly is the only contri- 
bution to S(x), and its matrix elements are renormali- 
zation group invariant. With Wilson fermions however, 
the absence of chiral symmetry implies that the trace 
anomaly acquires a linearly divergent contribution from 
the operator jpi/). Thus our matrix elements & g baro ' ) should 
be regarded as intermediate results. The coefficient of 
the counterterm, as well as its disconnected diagrams, 
will have to be computed before we can quote a physical 
value for b g in the pion. Not surprisingly, 6 g bare ^ shows a 
strong quark-mass dependence, since the missing discon- 
nected diagrams are suppressed by 1/m. We note that 
^(baic) ^ 0.9(1) at the largest mass is of the same order 
of magnitude as Ji's phenomenological estimate of b g in 
the proton [J, 0.85(5). 

Conclusion. — We have computed the glue momen- 
tum fraction (x) g in a pion of mass 0.6GeV < M v < 



1.06GeV using quenched lattice QCD simulations. We 
find 37(8) (12)% at f^j^g = 2GeV, a result compatible 
with phenomenological determinations m, [24| ■ 

Although it appears difficult to achieve precision at 
the percent level, the present method is applicable to full 
QCD with dynamical quarks. Presently the larger uncer- 
tainty comes from the normalization of the quark contri- 
bution to the renormalized (x) g , and could be reduced 
significantly by a one-loop calculation. 

We also evaluated the bare trace anomaly contribution 
to the pion's mass in the same framework. The countert- 
erm remains to be calculated, but it will ultimately be 
preferable to use chiral fermions to avoid mixing with the 
lower dimensional fermion operator. 

Finally, we remark that the freedom of choosing a 
numerically advantageous discretization of T M „ has not 
been fully exploited in previous lattice simulations. The 
improvement that was essential in the present computa- 
tion of the pion momentum fraction can be carried over to 
fully dynamical calculations and the exploration of other 
observables, such as the gluon contribution to the nucleon 
spin. It is also particularly promising for thermodynamic 
studies of pressure, energy density and transport coeffi- 
cients. 
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